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Femtosecond pulses in a dense two-level medium: Spectral transformations, transient 

processes, and collisional dynamics 

Denis V. Novitsk5l3 

B.I. Stepanov Institute of Physics, National Academy of Sciences of Belarus, 

Nezavisimosti Avenue 68, 220072 Minsk, Belarus 

Propagation of ultrashort optical pulses in a dense resonant medium is considered in the semi- 
classical limit. In our analysis, we place emphasis on several main points. First, we study trans- 
formations of spectra in the process of pulse propagation and interactions with another pulse. The 
second point involves the transient processes (including pulse compression) connected with self- 
induced transparency soliton formation inside the medium. Finally, the third aspect is the study of 
collisions of co- and counter-propagating pulses in the medium. In the last case, the investigation of 
symmetric and asymmetric collisions shows the possibility of effectively controlling the parameters 
of transmitted radiation. 

PACS numbers: 42.50. Md, 42.65. Sf, 42.65. Tg 



I. INTRODUCTION 

A two-level medium, i.e. a m.ediuin composed of a 
collection of two-level systems (atoms, molecules, etc.), 
is the basic quantum-mechanical model in the theory of 
light-matter interaction. It plays the same fundamental 
role as a harmonic oscillator in classical physics. There- 
fore, the comprehensive study of this model and its gen- 
eralizations is one of the most important problems in the- 
oretical optics. The two-level model has allowed the de- 
scription of a number of well-known effects IllH includ- 
ing the effect of self- induced transparency [J, |3| • Among 
the further generalizations of the model we single out the 
one which takes into account the near dipole-dipole inter- 
actions between the elements of the dense medium due 
to introduction into the equations of the so-called local- 
field correction [5|, |6[. In stationary regime, the presence 
of the local field results in the remarkable effect of intrin- 
sic optical bistability |7Hl0l|. 

In the regime of optical pulses, the investigation 
of local-field effects is linked with solitonic dynamics 
11, [ij], optical switching [ij, [ij], etc. It was shown 
14| that, in the dense two-level medium, the scale of 
population and field changes becomes so short that one 
cannot use the common slowly varying envelope approx- 
imation in space (SVEAS) to describe light pulse prop- 
agation. Therefore, one should use the full (not trun- 
cated) Maxwell wave equation. This is the case in the 
present paper where we consider one-dimensional prob- 
lem of ultrashort (femtosecond) light pulse propagation 
in a two-level medium. Two remarks should be made 
here. First, as discussed in our previous work [I^, the 
local-field effects depend strongly on the duration of the 
pulse. For the femtosecond pulses considered in this pa- 
per, these effects are negligible. However, the dynamics 
of the pulse cannot be fully described by the usual self- 
induced transparency (SIT) model. This fact leads us to 
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the second remark. In the dense two-level medium, cal- 
culations using the wave equation without SVEAS show 
the slow attenuation of SIT solitons as was predicted in 
Ref. [lal- This and similar other effects provide the rich 
dynamics of transient processes which are one of the main 
objects of this investigation. 

Another aim of our research is to study in detail the 
complex dynamics of pulse collisions in the two-level 
medium. Interaction of counter-propagating solitons was 
previously a research issue in both the cases of two-level 
[ITI [18| and three-level media j^lBi^l;] . In the appropri- 
ate section, we will return to some of this works and 
discuss the further progress achieved in this paper. Here, 
it is worth noting that we try to trace step-by-step the 
changes in parameters of the pulses obtained at relatively 
short distances due to the high density of the medium. 
In addition, collisions of solitons and attendant effects 
were theoretically and experimentally studied in different 
other systems such as optical fibers |23| , photorefractive 
crystals [2J, [23 , and photonic crystals [2^, Um ■ ^^ ^^^^ 
paper we consider temporal solitons, while the collisions 
between spatial ones were reported as well [22, [23 ■ 

The structure of the paper corresponds to the grad- 
ual transition from the relatively simple situations to the 
more complicated ones. Section II is devoted to the case 
of single pulse propagation in a dense two-level medium 
paying attention to the transformations of light spec- 
trum as it moves through the medium. In Section HI 
we consider interaction of co-propagating solitons, espe- 
cially those overtaking one another. Section IV is dedi- 
cated to the case of counter-propagating pulses colliding 
in the medium. In particular, we analyze the situation of 
the so-called symmetric collision there. Finally, in Sec- 
tion V we study the asymmetric collisions between the 
counter-propagating solitons. 



II. SINGLE PULSE PROPAGATION 

We start with the semiclassical Maxwell-Bloch system 
for population difference W, microscopic polarization R, 



and electric field amplitude J7' = fl/o; — {yujhbj)E (i.e. 
normalized Rabi frequency) [^, [2^ , 



—- = {"^W + iR{5 + tW) - 72i?, 
dr 

^ = 2i(r!'*ii' - R*9!) - 7i(P^ - 1), 
dr 



d^^' d^n 



,2o' 



9^2 
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where t = ujt and ^ = fcz are dimensionless argu- 
ments; S — Aid /uj is the normalized detuning of the 
field carrier (central) frequency w from atomic resonance; 
71 — {ujTi)~^ and 72 = (wT2)~^ are the rates of longi- 
tudinal and transverse relaxation, respectively; k = uj/c 
is the wavenumber, and c is the light speed in vacuum; 
e = ujl/oj = A-K fj? C / ihuj is the normalized Lorentz fre- 
quency. Here we assume that the background dielectric 
permittivity of the medium is unity (two-level atoms in 
vacuum). Note that in Eq. ([3]) we do not use slowly- 
varying envelope approximation (SVEA) which cannot 
hold true even for thin films of the medium [lj| . 

In this paper we consider propagation of ultrashort 
pulses with Gaussian shape f2 = r2pexp(— ^2/2^^) where 
tp is the pulse duration. Amplitude of pulses is mea- 
sured in the units of the characteristic Rabi frequency 
r^o = •\/27r/2ip, which corresponds to the so-called 2it- 
pulse [16|. For calculations the values Ti = 1 ns and 
T2 = 0.1 ns are taken, so that femtosecond pulses ap- 
pear to be in the regime of coherent interaction with the 
resonant medium. The spectra of pulses plotted in this 
paper are obtained as the absolute values of the Fourier 
transform of the corresponding field profiles. The spec- 
tra are normalized on the peak value of the incident 
pulse spectrum which is recognized as the unity. We use 
the following parameters of calculations which hold true 
throughout the paper if the other is not stated: 714 = 1, 
ujL = 10" s-\ 5 = 0, A = 0.5 fim, tp = 50 fs. 

In Fig. [1] we demonstrate the results of propaga- 
tion simulation for the pulse with high peak amplitude 
(Op = l.Sflo)- As it was shown in Ref. [la], an ultrashort 
pulse in dense resonant medium experiences compression. 
The reason for this is perhaps the effect of self-phase mod- 
ulation reported in Ref. [lj| . The effect of compression is 
characterized by a certain distance of optimal compres- 
sion: After this distance the pulse slowly attenuates due 
to dispersion. For the layer of approximately optimal 
thickness, the main part of the pulse is transmitted and 
compressed, while some part of its energy is absorbed 
and then reemitted [Fig. [Ua)]. The spectrum of the 
main part has the typical bell shape, while the reemitted 
light has a dip at the resonant wavelength [Fig. Hfd)]. 
This can be connected with effective reabsorption of low- 
intensity light at the resonant frequency. At the same 
time the main, high-intensity part of the pulse is almost 
not absorbed due to the self-induced transparency (SIT) 
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FIG. 1. (Color online) (a) Intensity profile of the pulse trans- 
mitted through the layer of thickness L = lOOA. (b) and 
(c) Spectra of transmitted and reflected light respectively. 
The spectra were calculated for the radiation appeared in the 
time interval from t = — 3ip (start of the incident pulse) to 
t — lOOtp. (d) Spectra of transmitted light calculated for the 
radiation appeared in the other time intervals (from t — — 3ip 
to i = 6tp, and from t — 6tp 
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FIG. 2. (a) and (b) Spectra of transmitted light in the cases 
of absence and presence of local-fleld correction, respectively. 
The spectra were calculated for the radiation appeared in 
the time interval from i — — 3tp (start of the incident pulse) 
to i = lOOip. (c) Difference between the spectra shown in 
plots (a) and (b). (d) Difference between intensity profiles 
of transmitted light in the cases of absence and presence of 
local- field correction, respectively. The amplitude of the pulse 
Qp — 1.5fio- 



mechanism. The resulting spectrum of the transmitted 
radiation is shown in Fig. [ijb): The two peaks are sit- 
uated at both sides of the central wavelength. The low- 
level reflected light has a dip in its spectrum as well [Fig. 
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FIG. 3. (Color online) (a) Intensity profiles of the light trans- 
mitted through the layer of thickness L = lOOA for the pulses 
with different peak amplitudes, (b), (c), (d) Spectra of trans- 
mitted light for the pulses with different peak amplitudes 
(yip = 0.5f2o, f^o and l.SJlo, respectively). The spectra were 
calculated for the radiation appeared in the time interval from 
t = —3tp (start of the incident pulse) to f = lOOtp. The 
dashed line corresponds to the spectrum of incident pulse. 
The local-field correction is taken into account. 
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FIG. 4. (Color online) (a) Intensity profiles of the light trans- 
mitted through the layer of thickness L = 500A for the pulses 
with different peak amplitudes, (b) Spectra of transmitted 
and reflected light for the pulse with the peak amplitude 
Qp — O.Sr^o- (c), (d) Spectra of transmitted light for the 
pulses with the peak amplitudes Jlp — flo and 1.5fio, respec- 
tively. The spectra were calculated for the radiation appeared 
in the time interval from t — —3tp (start of the incident pulse) 
to (c), (d) t = 200tp, (b) t = 600ip. The dashed line corre- 
sponds to the spectrum of incident pulse. 



Figure [T] was obtained in the absence of the local-field 
correction (LFC), that is without the term eW in Eq. 
dl]). To evaluate the influence of local field, we calculate 
the difference between the transmitted light spectra ob- 
tained with and without LFC [Fig. [H^c)]. It is seen that 
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FIG. 5. (Color online) Spectra of transmitted light for the 
pulse with the initial peak amplitude f7p = 1.5flo after prop- 
agation of (a) L = lOOOA, (b) L = 2000A. (c) Spectra of the 
main part of transmitted pulse and (d) spectra of the "tail" 
for L = lOOOA and 2000A. 



lO.On 







>• 6. 
4. 
2. 
0. 



(a) 



1 St pulse 



;'; incident 



28 
24 
20 
16 

12 



2nd pulse 



(b) 


1st pulse 


incident 


2nd pulse 




3rd pulse 


■'' 1 '■ u ^ ^ 


. A . 




0.48 



0.49 0.50 0.51 

Wavelength (urn) 



0.48 



0.49 0.50 0.51 

Wavelength (|j.m) 



0.52 



FIG. 6. (Color online) (a), (b) Intensity profiles of the light 
transmitted through the layer of thickness L = lOOOA for the 
pulses with peak amplitudes Qp — 2r2o and 3fio, respectively, 
(c), (d) Spectra of transmitted light corresponding to the up- 
per figures. 



LFC leads to the significant change in spectrum (more 
than 8% of the maximal value) , while the intensity pro- 
file remains almost unchanged [Fig. [5M)] in accordance 
with the results previously reported [l^- The difference 
between spectra is mainly due to the late, reemitted ra- 
diation which interacts with the medium for a time long 
enough to feel the local field. Nevertheless the overall 
form of the spectrum is the same as one can ascertain 
comparing Figs. [2][a) and (b). Therefore, hereinafter all 
the results are obtained taking into account LFC. 

Returning to Fig. [Jl we should once more emphasize 
the role of pulse "tail" in formation of the fine struc- 
ture of the complex spectrum of transmitted light. Note 
that the resulting spectrum is not a simple superposition 



of the two others [see Fig. [Ijd)] because of interference 
term which should be taken into account here. Figure 
|3] shows the spectra for the transmitted pulses with dif- 
ferent initial peak amplitude flp. The simplest case is 
observed for flp — flo, i.e. for the 27r-pulse. Such a 
pulse is only slightly compressed and can propagate for 
a long distance almost without change. Its spectrum is 
generally coincident with that of the incident pulse [Fig. 
ISlJc)]. The pulses with amplitudes larger and smaller 
than rig demonstrate significant pulse transformations. 
Low-intensity pulse (in our case flp = 0.5J7o) experiences 
strong absorption on the central frequency, so that the 
output radiation has no any regular profile [it is not seen 
in Fig. |3l[a) due to its proximity to zero]. The overall 
energy leaving the medium in this case is about half of 
the initial pulse energy (for the calculational conditions 
denoted in the caption of Fig. |3]). This output energy 
is mainly concentrated in the side bands of the spectrum 
as seen in Fig. ^h) . In the case of high-intensity pulse 
(Op = l.SOo), absorption is not significant and almost 
all energy transmits through the medium. However, its 
spectrum [Fig. ^d)] implies that some part of the ra- 
diation is transferred from the central wavelength to the 
side bands. 

The next step is to consider the transformations of 
spectra with distance traveled by the pulses in the two- 
level medium (we increase the distance by the factor of 5, 
i.e. to L = 500 A). We find that the spectrum of 27r-pulse 
gets more distorted as it propagates inside the medium 
for a longer distance [Fig. Hl^c)]. In Fig. Ht^b) one can see 
that the spectrum of low-intensity transmitted pulse has 
the same form as previously. In this figure the spectrum 
of reflected light is plotted as well. Reflected radiation is 
concentrated mainly very close to the central (resonant) 
frequency. In other words, the medium in this case can 
serve as a very-narrow-bandwidth filter. However, the 
efficiency of such a filter is not large because only a few 
percent of the incident light is reflected. 

The case Qp = 1.5^0 is of maximal interest to us due 
to the strongly pronounced role of the pulse "tail" . Com- 
parison of Figs. M^d) andSlJd) shows that the side bands 
move away from the resonant wavelength as the pulse 
propagates in the medium. This is also proved in Fig. 
[5ja) and (b). This transformation of spectrum is due 
to change in the spectrum of the "tail": The central 
dip gets wider [Fig. [IJd)]. Note that the main part of 
the pulse is characterized by the same spectrum [Fig. 
[SJc)]. It remains almost unchanged because the main 
part demonstrates invariant propagation in the regime 
of self-induced transparency (SIT-soliton). However, the 
main pulse very slowly attenuates though not so fast as 
was reported previously [16| . This is due to higher accu- 
racy of our calculations in the current investigation. 

If the pulse amplitude is high enough (Op > 20o), 
the splitting of the pulse into several solitons can be ob- 
served. The examples of spectral transformations con- 
nected with this splitting are shown in Fig. [S] Com- 
parison between the smooth spectra of separate solitons 



and the resulting jagged one implies that the phase rela- 
tions between frequency components plays an important 
role. This interference is especially characteristic for the 
pulse splitting into three components (Op — SOq). In 
this case [Fig. [DJd)], apart from the central band, the 
spectrum has several pronounced and wide side bands. 
It is also worth noting that the first soliton has always 
the widest spectrum while every next one is more narrow. 
This corresponds to the fact that the first pulse is more 
compressed than any consecutive one. 



III. CO-PROPAGATING PULSES 

Now let us consider the process of interaction between 
two pulses co-propagating in the dense two-level medium. 
Some preliminary results were published previously [16| . 
However, in this paper we present extended and more ac- 
curate calculational results. First, in Fig. [7l[a) we show 
how the intensity of the second pulse changes under the 
influence of the flrst one. The time interval between the 
peaks of two incident co-propagating pulses is taken to 
be 6ip. It is seen that the first pulse loosing part of its 
energy causes some increase in the intensity of the sec- 
ond one in comparison with the case of the single pulse 
(this corresponds to the dotted lines). This is especially 
evident in the case of Fig. [71[a4). Figures [TJb) demon- 
strate the transformations of spectra due to simultaneous 
compression and interactions between the pulses. Note 
that the spectrum of the pair O.SOq -t- IOq [Fig. [TJbl)] 
is similar to the spectrum of the single pulse with ampli- 
tude I.SOq [Fig. mjd)] and can be treated as composed 
from the spectra of pulses with O.50o and Oo [Fig. IH^b) 
and (c)]. Other spectra are characterized by appearance 
of strong side lobes and complex fine structure in the 
central part as seen in Fig. [7l^b2). 

One particular case (the pair 10o + 1.50o) is worth con- 
sidering in detail here. This is because the second pulse 
of this pair propagates in the medium faster than the first 
one as results from Fig. Hl^a). This means that, though 
the 27r-pulse does not infiuence the second pulse at short 
distances, their interaction should be getting stronger as 
they co-propagate in the medium. The reason for this 
distinction in pulse velocity can be easily understood on 
the basis of the well-known expression for the stationary 
pulse [i|. 
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where u is the group velocity of the pulse, c is the vacuum 
light speed, a ^ Ctp is the extinction coefficient, C is 
the concentration of two-level atoms, tp is the stationary 
pulse duration. Here we are interested only in relative 
speed of both pulses, so one can write 
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FIG. 7. (Color online) (a) Intensity profiles and (b) spectra of the light transmitted through the layer of thickness L — lOOA 
for the pairs of incident pulses. The numbers at (a) and (b) correspond to different peak amplitudes of incident pairs: (1) 
0.5Slo + f^o, (2) 0.5^0 + l-5ilo, (3) 1.5^0 + ^o, (4) l.Silo + l.bQo- The time interval between the maxima of incident pulses is 
6tp. The dotted line denotes (a) the transmitted single pulse with initial amplitude Qq (1) and l.Sr^o (2); (b) the spectra of 
initial pair of pulses. 
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FIG. 8. (a) Intensity profiles and (b) spectra of the light transmitted through the layers of different thicknesses. The calculations 
were performed for the pair of incident pulses fio + I-SOq. The time interval between the maxima of incident pulses is 6tp. 



or, taking into account that the peak intensity Ip ~ 1/^p, 
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This condition is perfectly satisfied as seen in Fig. HJa) 
(for the case Ci = C2). The concentration dependence 
also holds true as demonstrated by our careful exami- 
nations. This fact can be treated as one more proof of 
validity of our calculational scheme. 

Thus, we can study the process of collision of two 
pulses co-propagating with different velocities. This case 
is demonstrated in Fig. |8l It is seen that the second, 
more intensive pulse overtakes the first one and, finally, 
at the length approximately L = 250A they form a single 
pulse. This is indicated by the dramatic simplification 



of the spectrum which retains only one, central lobe. As 
the distance increases, this single pulse disintegrates, and 
the intensive pulse leaves the other behind. This fact 
corresponds to significant increase in complexity of the 
spectrum. The results of Fig. |8]also show that the inten- 
sity of the single pulse after coalescence is intermediate 
between the intensities of the pulses before the collision. 
The change in the peak intensities of the pulses as they 
co-propagate in the medium is depicted in Fig. |9l It is 
interesting to compare them with the behavior of inten- 
sity of a single pulse. It is seen that, after some transient 
process, the peak intensity of the single pulse remains al- 
most invariable. This corresponds to the regime of soli- 
tonic propagation with the envelope perfectly described 
by hyperbolic secant function (this means that 1.5f2o- 
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FIG. 9. (Color online) Dependence of peak intensity of the 
light transmitted through the layer on its thickness. The 
calculations were performed for the pair of incident pulses 
Qo + l.SSlo and for the single pulses Qq and l.Sr^o. 
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FIG. 10. Intensity profiles of the light transmitted through 
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identical envelope. The calculations were performed for the 
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pulse with initial area of Stt transforms into a common 
27r-soliton) . However, the intensity of the pulse is not 
strictly constant, it slightly fluctuates; so we cannot talk 
about strict, mathematical solitons but rather about re- 
alistic, physical ones. In addition, the amplitude of these 
fluctuations tends to decrease with the propagation dis- 
tance. However, at larger distances the pulse slowly at- 
tenuates. For example, the intensity of the pulse with 
initial amplitude l.SrJo drops from approximately A.BHq 
(at L = lOOOA) to 3.75f]^ (at L = lOOOOA). 

The next question is connected with the nature of the 
transient process: At first, the pulse is strongly com- 
pressed (on the distance of about lOOA), and only then 
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FIG. 11. (Color online) Spectrum of light transmitted 
through the layer of thickness L — 1000 A. The calculations 
were performed for the pair of counter-propagating incident 
pulses with amplitude fio. 
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FIG. 12. (Color online) Dependence of peak intensity of the 
light transmitted through the layer on its thickness. The cal- 
culations correspond to the case represented in Fig. [lOl 



its envelope transforms to the solitonic one (hyperbolic 
secant) which is connected with some decrease in the 
pulse intensity. The initial compression of the pulse is 
due to the self-phase modulation effect and is strongly 
dependent on the form of the incident pulse. This com- 
pression can be observed even for the pulses with initial 
envelope described by hyperbolic secant. It is known 
that, at the same area (say, 27r), the sech-pulse is wider 
than the Gaussian one, so that if the soliton is formed 
from the initial pulse, it would obtain greater duration. 
But it is not the case: the soliton forms from the com- 
pressed pulse. This compression can be characterized by 
the length value, namely, the so-called distance of opti- 
mal compression. 

Returning to the collision of two pulses, we see (Fig. ^ 
that they strongly influence one another, but, after the 
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FIG. 13. Distribution of (a) population difference and (b) light intensity inside the layer L — lOOOA = 500 /im at different time 
points. The calculations were performed for the pair of counter-propagating incident pulses with amplitude r^o- 



change of the propagation order, the distance between 
them increases and the sohtons are restored. One can say, 
that this figure demonstrates the interaction between the 
already formed or just forming solitons which are colhsion 
resistant. 



IV. COUNTER-PROPAGATING PULSES 

From the previous section we know that the collision 
of co-propagating pulses does not prevent formation of 
stationary pulses with the same characteristics as in the 
single pulse case. The collision of counter-propagating 
pulses has fundamentally different results. The reasons 
were studied previously in Refs. [13, [l8|. From math- 
ematical point of view, two counter-propagating pulses 
cannot provide a stationary (solitonic) solution. Phys- 
ically, this means that the absorption of some part of 
energy of the pulses occurs in the region of collision. Fur- 
ther we consider the dynamics of such collisions in detail. 

We start from the collision of two Gaussian pulses with 
the amplitude f2o entering the medium at the same in- 
stant of time. In such symmetric geometry the light ap- 
pearing from both ends of the layer (we can conditionally 
call it transmitted and reflected light) will be described 
by the same envelope. These envelopes are shown in Fig. 
[TUlfor different thicknesses of the layer. It is seen that the 
result of the collision strongly depends on the thickness 
L and, as L increases, the long precursor and tail be- 
come apparent (in Ref. [lM| only the tail was observed). 
The spectra of the precursor and the tail are similar (Fig. 
Ilip : they have a dip in the region of resonant wavelength, 
while radiation at this central frequency is strongly ab- 
sorbed. 

To study this thickness dependence of the collision, we 
simulate the process of interaction between the pulses in 
a wide range of thicknesses L. The results of these cal- 



culations for the peak intensity of the transmitted pulse 
are demonstrated in Fig. [T^l It is seen that at relatively 
small length there are pronounced features which have a 
structure of something like resonances. For larger thick- 
nesses the amplitudes of these resonances diminish, so 
that the peak intensity of the pulse tends to the approx- 
imately constant value. This implies that the resonances 
observed are connected with some transient process. This 
process can be traced using distribution of population 
difference and light intensity inside the medium at dif- 
ferent instants of time (Fig. [T3l) . The thickness of the 
layer (L = lOOOA) is large enough for 27r-solitons to be 
formed by the time of collision. This stationary pulse is 
characterized by the typical population difference profile 
(drop to full inversion and following rise to the ground 
state). One can easily ascertain that, at the location of 
collision, a strong absorption really occurs, so that some 
part of light energy remains inside the medium even af- 
ter the pulses have gone away from each other. How- 
ever, after collision, pulse propagation is not simple: Af- 
ter reaching minimum at i = 32tp, the pulse intensity 
grows {t = 34ip), then goes down again {t = 36tp), and 
so on. This process is obviously connected with fluctu- 
ations depicted in Fig. [5] and results in precursor and 
tail formation. It is reasonable to suggest that the res- 
onances of Fig. [12] are due to pulse exiting in different 
moments of this transient process. But if the thickness of 
the layer is large, these fluctuations become smaller and 
smaller, so that at the output we have almost stationary 
pulse with the area 27r again and somewhat decreased 
intensity [t — 65tp). 

The next important feature of pulse collisions is their 
dependence on the duration (or, equivalently, peak inten- 
sity) of the initial pulses. This dependence is shown in 
Fig. [Mia) and is in accordance with the results of Ref. 
[rn |: For short and high-intensity J7o-pulses, the influence 
of collision on pulse propagation is not significant. On the 
contrary, long and low-intensity pulses strongly interact 
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with each other and, if the duration is large enough, the 
main pulse (soliton) disappears leaving only a precursor 
and a tail [see Fig. IWb) for the pulses with tp = 100 fs]. 
Note that the area of the main pulse (27r) is conserved 
even in the case of dramatic decrease of its peak inten- 
sity. For example, the pulse corresponding to tp = 70 fs 
in Fig. [HTb) can be easily fitted by a hyperbolic secant 
function. This means that such low-intensity solitons un- 
dergo substantial broadening. 



V. ASYMMETRIC COLLISIONS 

In previous section we considered the process of colli- 
sion of the two identical counter-propagating pulses en- 
tering the medium at the same instant of time. Such 
geometry of the problem can be referred to as the sym- 
metric one. Further we discuss the case of asymmetric 
collisions when the incidence time or intensity is different 
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for both initial pulses. 

Figure [TS] shows the results of collision of the pulse 
with the amplitude J7o with the identical pulse which en- 
ters the medium at different time. Negative values of 
the incidence time difference At means that the counter- 
propagating pulse is launched earlier than those propa- 
gating in the forward direction {At = corresponds to 
the symmetric situation considered in the previous sec- 
tion). In the case of At < 0, the collision occurs closer 
to the input of the forward-propagating pulse, so that 
it has much time to form a stationary pulse after the 



(a) 



0.5 
0.0 


^ 


t=5t^ 


N 


0.5 








1.0 









f=24f 




100 200 300 400 500 



1 00 200 300 400 500 

Distance (iim) 



1 00 200 300 400 500 



100 200 300 400 500 

Distance (nm) 



I 200 300 400 500 



FIG. 17. Distribution of (a) population difference and (b) light intensity inside the layer L = lOOOA at different time points. 
The calculations were performed for the pair of counter-propagating incident pulses: one with amplitude Q,o and the other with 
1.5^0- 



collision. Therefore, the fluctuations of the peak inten- 
sity [Fig. llSf a)] are less pronounced in this case than 
for At > 0. On both edges of this dependence (that 
is, for large absolute values of At) one obtains the triv- 
ial case of the single pulse propagation. The transition 
to this extreme regime has jump-like character on both 
sides of the peak intensity dependence. This symmetry 
does not take place for the time of pulse transmission 
through the layer [Fig. [TSl b')]. For large negative At 
the abrupt increase of time propagation occurs due to 
formation of the low-intensity (and, hence, slow) soliton 
near the very input. As the counter-propagating pulse 
enters the medium later and later, the fast initial pulse 
has enough time to pass a large fraction of distance be- 
fore the collision happens. As a result, the propagation 
time decrease gradually as seen in Fig. [T5lb). Thus, 
the coUisional dynamics discussed imply the possibility 
of controlling such pulse parameters as its peak intensity 
and transmission time by proper choice of incidence time 
of the counter-propagating pulse. 

Another approach to controlling the parameters of 
the transmitted pulse is to use the counter-propagating 
pulses of different intensities. In Fig. [TBI one can see the 
change in peak intensity of ilo-pulse as a function of the 
amplitude of the counter-propagating (controlling) one. 
As this amplitude grows, the peak intensity of the trans- 
mitted forward-propagating pulse decreases as well as its 
propagation velocity. Finally, in the region of amplitudes 
near LSf^Oj the transmitted solitonic pulse is absent per 
se. We have only the precursor and the tail. For larger 
amplitudes of controlling pulse, the transmitted one ap- 
pears again, of course, with larger retardation. However, 
coUisional dynamics in this case are even more complex 
due to the effects of pulse splitting. Therefore, we will 
not consider the details of these multi-pulse interactions 
in the present investigation. 

Instead, we turn to the case of complete soliton disap- 



pearance at the output of the medium when the ampli- 
tude of the controlling pulse equals 1.5f2o- The dynamics 
of level population difference and light intensity inside 
the medium depicted in Fig. [T7] imply that the primary 
importance in this effect is likely to be connected with 
the residual excitation of the medium after passing the 
controlling pulse. Since its initial area is Stt (amplitude 
1.5rJo)j before the 27r-soliton formes, some part of the 
energy is absorbed giving rise to the pronounced tail. It 
is seen that the collision leads to the significant loss of 
energy of the forward-propagating pulse while the con- 
trolling (high-intensity) pulse remains almost unchanged. 
Moving farther, the low-intensive pulse meets the tail 
of the counter-propagating one and the residual excita- 
tion and, finally, undergoes ultimate absorption in the 
medium. This absorbed energy gives fluorescent radia- 
tion in the long run as was stated in Ref. |18| . 



VI. CONCLUSION 

In this paper, we have considered the ultrashort (fem- 
tosecond) pulse interaction with a dense collection of 
two-level atoms. In other words, the light-medium in- 
teraction was investigated in the coherent regime when 
incoherent phenomenological relaxation can be neglected. 
Our study was conducted in semiclassical approximation, 
which is valid as we are not interested in the detailed de- 
scription of the processes of spontaneous radiation [l|. 
In addition, semiclassical approach allows us to correctly 
take into account the local field correction [6] . Our calcu- 
lations were performed for two-level atoms located in vac- 
uum which corresponds to dense gaseous media. In order 
to consider solid-state systems (two-level atoms embed- 
ded in a dielectric, for example, quantum dots in semi- 
conductor or glassy environment), one has to take into 
account the so-called local-field enhancement factor, as 
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it was done in stationary regime [29|, |30|. However, if the 
background dielectric is not absorptive, behavior of op- 
tical pulses reported in this paper is expected to be still 
valid for this more complex case (with the corresponding 
renormalization of the electric field) . 

In conclusion, we should discuss the prospects of pulse 
collision dynamics studied in the paper. Our accu- 
rate calculations show the self-healing effect for the co- 
propagating solitons but this is not the case for the 
counter-propagating pulses. Together with the transient 
processes of soliton formation, the strong interaction of 
the counter-propagating pulses provides for the rich dy- 
namics of population difference and light intensity. The 
proper choice of the parameters of the colliding pulses al- 
lows one to effectively control the characteristics of light 



at the output of the medium. For example, one can ob- 
tain almost entire absorption of the soliton whose en- 
ergy can be stored inside the medium, at least, for the 
relaxation time. The conditions of controlling pulse in- 
tensity by choosing the intensity or incidence time of the 
counter-propagating one can be considered as a basis of 
the peculiar logic gates or other elements for which the 
possibility of realization is still to be studied. 
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